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^ We present parsec-scale resolution observations of the atomic and molecular ISM in 

I 1 two Galactic supershells, GSH 287+04-17 and GSH 277+00+36. Hi synthesis images 

^ from the Australia Telescope Compact Array are combined with ^^CO(J=1-0) data 

p> from the NANTEN telescope to reveal substantial quantities of molecular gas closely 

^ associated with both shells. These data allow us to confirm an enhanced level of molec- 

ularization over the volumes of both objects, providing the first direct observational 
evidence of increased molecular cloud production due to the influence of supershells. 
We find that the atomic shell walls are dominated by cold gas with estimated temper- 
O atures and densities of T ^ 100 K and no ^ 10 cm~^. Locally, the shells show rich 

substructure in both tracers, with molecular gas seen elongated along the inner edges 
.Jph of the atomic walls, embedded within Hi filaments and clouds, or taking the form of 

^ small CO clouds at the tips of tapering atomic 'fingers'. We discuss these structures in 

the context of different formation scenarios, suggesting that molecular gas embedded 
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within shell walls is well explained by in-situ formation from the swept up medium, 
whereas CO seen at the ends of fingers of Hi may trace remnants of molecular clouds 
that pre-date the shells. A preliminary assessment of star formation activity within the 
shells confirms ongoing star formation in the molecular gas of both GSH 287+04-17 
and GSH 277+00+36. 

Subject headings: ISM: atoms, ISM: bubbles, ISM: clouds, ISM: evolution, ISM: molecules, 
radio lines: ISM 



Introduction 



The neutral interstellar medium (ISM) of the Galactic disk is riddled with loops, shells and 



cavities (e.g. |Heiles||1979| [McClure- Griffiths et ai:]|2002t [Ehlerova fc Palous||2005| ). The largest of 
these structures - Hi supershells - may reach hundreds of parsecs in diameter, with formation 
energies as high as ^ 10^^ ergs. Such objects strongly influence the structure and evolution of the 
Disk ISM. 

The dominant paradigm holds that supershells are formed through the cumulative action of 
multiple stellar winds and supernovae, which blow hot, overpressurized bubbles, and sweep up the 
surrounding medium into cool, dense shells ( Bruhweiler et al.|1980 ; Tomisaka et al.|1981 ; McCray &; 
Kafatos|1987 ). The accumulation of the ISM in such superstructures is one means of generating the 
high densities and column densities required for the production of molecular gas, and supershells 



have long been suggested as drivers of molecular cloud formation (e.g. McCray & Kafatos 1987 



Mashchenko & Silich|1994t Fukui et al.|1999; 


Hartmann et al.|2001|). However, despite a substantial 


body of theoretical work (e.g. 


Koyama &; Inutsuka 


200C 


); Bergin et al. 


2004; 


Heitsch &; Hartmann 


2008; 


Inoue &; Inutsuka||2009), and an growing list of articles documenting supershell-associated 


molecular gas (e.g. 


Handa et al. 1986 


; Jung et al. 1996; 


Fukui et al. 1999; 


Matsunaga et al. 2001; 


Yamaguchi et al.|2001; Dawson et al.|2008a 


), conclusive observational evidence of this phase change 



occurring in the walls of shells has not yet been found. The degree to which large-scale stellar 
feedback is driving the production of molecular clouds remains unconstrained. 

Conversely, in a pre-structured and highly inhomogeneous ISM, an expanding supershell will 
undergo encounters with pre-existing dense clouds, which may compress, disrupt, fragment, or 



even completely destroy them (e.g. Klein et al. 1994; Foster &; Boss 1996; Mellema et al. 2002). 
The relationship between supershells and the molecular ISM should therefore be governed by the 
interplay between this potentially destructive process and the creation of new molecular material. 

This paper presents parsec-scale resolution observations of atomic hydrogen and ^^CO(J=1-0) 
in two Galactic shells, GSH 287+04-17 and GSH 277+00+36, with the aim of investigating the 
role played by supershells in the evolution of the molecular ISM. These are some of the highest 
resolution images of any supershell, and the first time a dedicated comparison between the atomic 
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and molecular material in shell walls has been performed. 



In the remainder of this introduction we briefly summarize the basic properties of the two shells, 
before moving on to describe the observational and data reduction techniques in ^ Sections |3.1 



and |3.2| discuss the distribution and morphology of the two tracers, and examine the properties of 



the two phases quantitatively via Gaussian fitting to the shell spectra. Section [3^ investigates the 
degree of molecularization in the shell volumes, finding evidence of enhanced molecular fractions 



in both objects, and §3.4| gives their total Hi and H2 masses. In §4.1| and §4.2| we discuss Hl-CO 
structures in the shell walls in the context of different formation scenarios. Section 14.31 takes a 
preliminary look at star formation activity in the shell molecular clouds, and our conclusions are 
finally summarized in ^ 



1.1. GSH 287+04-17: Basic Properties 



GSH 287+04-17, also known as the 'Carina Flare' supershell, is a medium-sized, gently ex- 
panding Galactic chimney, located in the Sagittarius-Carina Arm at a distance of 2.6 ± 0.4 kpc 
and Galactocentric radius of ^ 8 kpc. It was originally discovered in ^^CO(J=1-0) by Fukui et al. 



(1999), who reported a scattering of molecular clouds extending in a wide swath above the Galactic 



Plane. These clouds showed telltale signs of global expansion, leading the authors to correctly 



identify them as parts of an expanding superstructure. Dawson et al. (2008b) later confirmed the 



existence of the counterpart atomic shell using low resolution ^ 16' Hi data from the Southern 



Galactic Plane Survey (SGPS; McClure-Griffiths et al. 2005), revealing its global properties and 
large-scale morphology for the first time. The main body of the shell measures ^ 230 x 360 pc, and 
it has broken out of the disk at a height of z ^ 280 pc, with a associated high-latitude emission 
seen up to heights of ^ 450 pc above the midplane. The molecular clouds form co-moving parts 
of the shell, which is estimated to contain total Hi and H2 masses of Mhi - 7 ± 3 X 10^ Mq and 
^ 2.0 + 0.6 X 10^ M0, respectively. The expansion velocity of the shell is ^ 10 km s~^, and its 
age and formation energy are estimated to be ^ 1 x 10^ yr and between 0.5 — 1 x 10^^ ergs, based 
on comparisons with analytical and numerical models. 



1.2. GSH 277+00+36: Basic Properties 

GSH 277+00+36 is a large outer Galaxy chimney, located at the edge of the Sagittarius-Carina 
spiral arm at a kinematic distance of 6.5 ± 0.9 kpc and Galactocentric radius of ^ 10 kpc. It was 



originally discovered in the low resolution portion of the SGPS (McClure-Griffiths et al. 2000), 
where it is seen as a prominent void in the bright Hi emission of the Plane. With a main body 
^ 610 pc in diameter and chimney extensions that reach more than 1 kpc above the midplane, this 
large and evolved supershell has swept up an estimated atomic mass of Mhi ~ 3 ± 1 X 10*^ M©, 
and is expanding with a velocity of ~ 20 km s~^. The age and formation energy of the shell are 
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estimated to be between 1 — 2 x 10^ yr and 0.9 — 2.4 x 10^^ ergs, where the higher values in these 
ranges are taken directly from McClure-Griffiths et al. (2000) and the lower values are derived using 
the same analytical formula as Dawson et al.| ( |2008b ) in order to better illustrate the difference 
between the two shells. High-resolution 3') Hi observations by McClure-Griffiths et al. (2003) 
revealed narrow, well-defined walls, sharply delineated along their inner edges, with a great deal of 
complex substructure including knots, filaments and 'drips'. No molecular observations have been 
previously reported for this shell. 



2. Observations 
2.1. HI Observations and Reduction 



This paper makes use of new Australia Telescope Compact Array (ATCA) mosaicing obser- 
vations of GSH 287+04-17 and previous observations of GSH 277+00+36 from [McClure-Griffiths 



et al. (2003). Here, the GSH 277+00+36 mosaicing observations are reprocessed completely from 



the original raw data files to ensure consistency in the data reduction procedures for both shells. 

GSH 287+04-17 was observed over two sessions in June and November of 2007 with the EW352 
and EW367 array configurations. These configurations provided nearly uniform baseline coverage 
from 31 to 367 m at intervals of approximately 15 m. The observations covered an 81 deg^ region 
between 283 < / < 292° and 1.5 < 6 < 10.5°, in 653 pointings arranged in a hexagonal pattern at a 
spacing of 23^ Although this value is slightly less than Nyquist spacing for the ATCA 34' primary 
beam, the sensitivity function only varies by ^ 2.5% between pointings, and the smaller number 
of pointings allows more snapshots per pointing to improve u-v coverage. A total of fifteen 40 s 
integrations were made for each pointing at widely spaced hour angles, and the standard sources 
PKS 1934-638 and PKS 1039-47 were observed for bandpass and absolute flux-density calibration, 
and amplitude and phase calibration respectively. The correlator configuration provided a velocity 
channel width of 0.41 km s~^ at 1.420 GHz. 

GSH 277+00+36 was observed in three sessions between October 2001 and March 2002. The 
observing strategy for GSH 277+00+36 differed from that of GSH 287+04-17 in the fohowing 
points: 1. A total of 841 pointings were made covering the regions 272 < / < 284°, 1.0 < 6 < 5.5° 
and 272 < / < 284°, —7.0 < 6 < — 1.0°; 2. each pointing was observed in eight 60 s snapshots; 
3. the standard sources used for amphtude and phase calibration were PKS B1039-47 and PKS 
B0843-54; 4. the correlator configuration provided a velocity resolution of 0.82 km s~^ at 1.420 
GHz. 

For both supershells, coverage of the Galactic Plane region {\h\ < 1.0°) was provided by ATCA 



data taken during the Southern Galactic Plane Survey (SGPS, McClure-Griffiths et al. 2005). 
Pointings were spaced 19' apart in a hexagonal configuration and a total of forty 30 s integrations 
were made towards each point. The array configurations used were 210, 375, 750A, 750B, 750C, and 
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750D, with PKS B1934-638 observed for bandpass and absolute flux calibration and PKS 0823-500 



observed for phase and amplitude calibration. Further details may be found in McClure-Griffiths 



et al. (2005). 



Calibration and imaging were performed using standard routines from the MIRIAD software 
package (|Sault fc Killeen|[2009[). Calibration, baseline subtraction and Doppler correction were 



carried out in the u-v domain. The individual pointings were then linearly combined and imaged 
using the MIRIAD task INVERT - a standard grid and fast Fourier transform technique - with 
superuniform weighting. Deconvolution was performed using the maximum entropy-based decon- 



volution algorithm MOSMEM (Sault et al. 1996), and the images were restored using the task 



RESTOR. For both shells the dimensions and position angle of the dirty beam varied somewhat 
with the pointing. For this reason we chose to restore with a circular gaussian beam whose FWHM 
matched the largest values measured for the major axes of the dirty beams. This corresponds to 
2.5' for GSH 287+04-17 and 3.0' for GSH 277+00+36. 

The synthesis images for both supershells were then combined with Parkes telescope single 
dish data from the Galactic All-Sky Survey (GASS, [McClure-Griffiths et al]|2009t [Kalberla et al. 
2010). GASS has an effective angular resolution of ~ 16', an effective velocity resolution of 1.0 
km s~^, and a brightness temperature sensitivity of ^ 60 mK. The GASS data were re-gridded to 
match the ATCA data, and the two were then linearly combined in the Fourier domain using the 
MIRIAD task IMMERGE. 

The final synthesized beam sizes are 2.5' and 3.0' for GSH 287+04-17 and GSH 277+00+36 
respectively, and the velocity resolution of both datasets is 1.0 km s~^ (although the channel width 
of the gridded data is 0.82 km s~^). The r.m.s. noise in a single 0.82 km s"-*^ channel is ^ 1 K for 
both objects. 



2.2. ^2CO(J=1-0) Observations and Reduction 



■'^^CO(J=1-0) observations of the GSH 287+04-17 region were made in two sessions from April 
to November 1998 and November 2001 to May 2002, using the NANTEN 4m telescope. Partial or 



complete datasets have previously been published in several studies (Fukui et al.||1999; Matsunaga 



2002; Dawson et al. 2008b), and most recently have been presented in catalogue form by Dawson 



et al. (2008a). The telescope half power beam width at 115 GHz is ^ 2.6'. All observations were 



made by position switching, with pointing centers arranged in a square grid at spacings of 2' or 4'. 
The vast majority of supershell-associated CO emission 99% of the total integrated intensity) is 
covered at 2' spacing. Further details of the observing strategy and region coverage may be found 
in 



Dawson et al. (2008a). 



Intensity calibration was made by the chopper wheel method, and Orion KL (q^bi950 — 
5^32^47.5^, 5^1950 = — 5°24'21") was observed as a standard calibrator source. The 2048 channel 
acousto-optical spectrograph provided a velocity coverage and resolution of 100 km s~^ and 0.1 
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km s ^, respectively. Typical system noise temperatures at 115 GHz were ^250 K, resulting in an 
r.m.s. noise of ^ 0.5 K per channel. 

The GSH 277+00+36 region was observed during the late 1990s as part of a pilot run for 



the NANTEN Galactic Plane Survey (Mizuno fc Fukui 2004). The data cover a Galactic latitude 



range of |6| < 5° between 272 </ < 280° and -4 < 6 < +3° between 280 </ < 284°, with 
velocity coverages of —100 < visr < 100 km s~^ and —50 < visr < 50 km s~^ for these two regions 
respectively. The observations were made by position switching at a grid spacing of 4', at a velocity 
resolution of 0.65 km s~^. The final data cube is gridded to a resolution of 1.0 km s~^. The r.m.s. 
noise per channel is typically ^ 0.35 K. 

The differing data quality in the two supershell regions necessitates further brief comment. In 



GSH 287+04-17 the detection limit in Leo is estimated to be ^ 3.5 pc^ K km s ^ (Dawson et al 



2008a), where Lco is the sum of the velocity-integrated intensities over the entire projected area of a 
cloud. For GSH 277+00+36, the lower sensitivity, wider grid spacing and larger distance, combined 
with non-Nyquist noise arising from poor baseline stability, conspire to produce a detection limit 
approximately an order of magnitude higher, which affects our ability to reliably detect small or 
low- brightness clouds. 



Section |3.3| briefly makes use of NANTEN Galactic Plane Survey datacubes covering regions 
immediately adjacent to the two shells. This additional data covers the regions 284 < / < 293°, 
-4 < 6 < +2° for GSH 287+04-17, and 270 </ < 272°, -5 < 6 < +4° for GSH 277+00+36. 
The former is a composite of observations made at grid spacings of either 2' or 4', at a velocity 
resolution of 0.1 km s~^, with r.m.s. noise levels per channel ranging between ^ 0.3 and ^ 0.9 K. 
The latter was made at a grid spacing of 8^ and has a velocity resolution of 1.0 km s~^, with an 
r.m.s. noise of ~ 0.2 K per channel. 



3. Results 

3.1. Observational Characteristics of Hi and CO in the Two Shells 

3.1.1. GSH 287+04-17 

Figure [l] shows velocity channel maps of the Hi and CO(J=1-0) emission in the GSH 287+04- 
17 region, both at a spatial resolution of ^ 2 pc. Each image displays an average over three Hi 
velocity channels, corresponding to a width of 2.46 km s~^. These new high-resolution Hi data 
reveal a wealth of clumpy, knotted and filamentary substructure in the atomic shell walls, on size- 
scales ranging down to the resolution limit. The global morphology of the supershell has been 



described in detail by Dawson et al. (2008b) and we do not duplicate their discussion here. Instead 
we focus on the information this matched-resolution Hi and CO data provides on the relationship 
between the atomic and molecular medium on local scales. 
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At velocities below ^ —20 km (panels a to / in figure [T]) velocity crowding is minimal, 
and emission belonging to the supershell is clear and unambiguous. The atomic and molecular 
ISM form striking configurations. These include cases in which molecular gas is elongated along 
the inside of the atomic shell wall, neatly embedded within Hi filaments or extended clouds, and 
numerous cases of small CO clouds at the tips of tapering Hi features, whose morphology suggests 
sculpting of the atomic medium by the shell energy source. 

Figure [2] shows a closeup of some of the most prominent CO-Hl features in the datacube. Panel 
al shows a section of the approaching wall; hereon referred to as the 'approaching limb complex'. 
Here, tapered fingers of Hi gas point radially inwards towards the shell interior, some tipped with 
clumps of CO-bright molecular gas. The shell walls between these features form arch-like structures, 
with a typical separation of ^ 20—40 pc. The CO clouds are small, with estimated molecular masses 
of 70 < < 600 Mq ( [Dawson et al.|2008aD . In contrast, the CO-Hi complex at (/, h) ^ (286°, 5°) 
shows a bright, elongated molecular cloud apparently embedded within the curved Hi wall. Here, 
the CO and Hi distributions are tightly correlated, to the extent that thin, twisted filaments within 
the wall are in places neatly traced in both transitions (shown in panel a2 of the figure). The mass 
of the largest CO cloud within this complex is estimated to be M^^ ^ 5500 M©, and it contains 
an active star forming region at (/,6) ^ (285.9°, 4.5°) ( [Dawson et al.|[2007( see also ^4.3). 

Panel b of the same figure shows emission originating from the high latitude regions of the 
supershell. The prominent group of CO clouds in this region are located z ^ 400 pc above the 
Galactic midplane. This is ^ 8 times higher than the molecular gas scale height at the Galacto- 
centric radius of the shell (a^ ^ 50 pc at i? ^ 8 kpc; Clemens et al. 1988; Malhotra 1994) and 
illustrates the likely importance of supershells in providing dense material to the Disk-Halo inter- 
face. This region, hereafter referred to as the 'high latitude complex', is dominated by a large, 
bright V-shaped atomic structure extending over 200 pc in length, whose foremost regions con- 
tain an estimated ^ 1.1 x 10^ Mq of molecular hydrogen. The leading CO cloud in this group is 
located at the very tip the atomic feature, offset from the bulk of the Hi emission. In contrast, 
the remainder of the molecular gas in the region is generally well associated with bright knots and 
filaments in the atomic medium. To the left of this complex, a single isolated molecular cloud at 
(/, b) ^ (290.3°, 7.6°) is seen at the very end of a long, filamentary arc of Hi. 

At velocities greater than —20 km s~^ velocity crowding begins to become problematic, and 
it is not always possible to disentangle shell emission from the bright and complex Hi background. 
Nevertheless, some associated structures are evident, and show similar characteristics to the features 
described above. Panel c shows a complex identified as part of the shell's receding limb, which 
contains several of the brightest molecular clouds in the sample. These clouds occupy positions 
along the lower edges of the Hi feature, on the side facing the Galactic plane, with a single cloud 
at (/,6) ^ (287.1°, 2.4°) offset towards the tip of a protrusion of atomic gas. The largest have 
estimated H2 masses of between several hundred and ^ 2000 Mq, and show signs of active star 
formation (see ^4.3). 
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3.1.2. GSH 277+00+36 

Our CO(J=1-0) data newly reveal a population of molecular clouds scattered throughout the 
walls of GSH 277+00+36. The global properties and detailed morphology of the atomic shell have 



been described previously by McClure- Griffiths et al. (2000, 2003), the latter using the same high- 



resolution Hi observations as the present work. However, the comparison with molecular tracers 
is entirely new, and sheds the first light on the hitherto overlooked molecular component of this 
evolved outer Galaxy shell. All H2 masses quoted in this section are derived using the X-factor as 
described in §3.2.2 



Figure [3] shows velocity channel maps of the Hi and CO(J=1-0) emission in GSH 277+00+36, 
at spatial resolutions of ^ 6 pc and ^ 7 pc, respectively. Each image displays an average over three 
Hi velocity channels, corresponding to a width of 2.46 km s~^. Molecular clouds are scattered 
throughout the main walls, forming co-moving parts of the Hi shell, with the clearest associations 
seen around the systemic velocity (panels e to j). Despite differences in age, distance and evo- 
lutionary stage, features in the shell walls show notable similarities to GSH 287+04-17, with CO 
clouds found embedded within the main walls, offset towards the tips of Hi features, and elevated 
to high altitudes above the midplane. 

Figure |4j panel al shows a closeup of a section of 'scalloped' Hi wall, in which the atomic walls 
form arcs and 'fingers' similar to those seen in the approaching limb of GSH 287+04-17 (panel al 
of figure [2]). A molecular cloud of M^^ ^ 1000 Mq is offset towards the tip of one of these atomic 
protrusions. To the right of this feature CO emission is detected in a single position emanating from 
the very tip of a smaller, narrower 'drip' of Hi protruding from the shell wall (panel a2). This small, 
unresolved CO cloud is on the very limit of detect ability for the dataset, permitting the possibility 
that other similar atomic 'drips' on the inner walls of the shell could harbor as- yet-undetected 
molecular components. This region is hereon referred to as the 'molecular drip' region. 

Panel b of the same figure shows a large, bright molecular cloud complex embedded within the 
main shell walls. This prominent group of CO clouds spans ^ 100 pc, with an estimated molecular 
mass of ^ 5 x 10^ Mq. Its rightmost member contains the giant Hll region RCW 42. 

Panel c shows a large (^ 5000 Mq) molecular cloud associated with the walls of one of the 
shell's chimney conduits. Located z ^ 450 pc below the Galactic midplane, its altitude is ^ 6 
times the molecular gas scale height in the outer Galaxy (cr^ ^ 80 pc; Heyer et al.||1998 ), providing 



another example of a massive high-latitude molecular cloud associated with a known supershell. 



3.1.3. Morphology Classes: Embedded and Offset CO Clouds 

The discussion above refers to 'embedded' CO clouds, that appear to be well embedded in parts 
of the atomic shell, and 'offset' CO clouds, that are found at the tips of Hi features, offset from 
the bulk of the atomic material. In the analysis to follow we find it useful to explicitly distinguish 
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between these two loose morphological categories. Those clouds treated as 'embedded' and 'offset' 
are marked in figures [2] and [4j 



3.2. Quantitative Properties from Gaussian Decomposition 

Gaussian decomposition of Hi spectra enables the separation shell emission from the multi- 
component background, and hence the extraction of accurate line profiles from which to obtain 
quantitative information on the atomic shell ISM. Equally importantly, it enables meaningful com- 
parisons with the ^^CO(J=1-0) line, allowing us to examine how Hi properties may vary according 
to the presence and characteristics of molecular gas. 

Here Gaussian fitting is performed for a sub-sample of spectra from the shells. We select 
the most prominent Hl-CO features in the approaching limb and high latitude complexes in GSH 
287+04-17 (panels a and b in Figure [2]), and the GSH 277+00+36 molecular 'drip' region (panel a 
of Figure |4|. Shell emission from these regions is generally distinct and unconfused, making them 
well suited to this quantitative analysis. The sample of observed spectra in the GSH 277+00+36 
region is small, and is included mainly to check for consistency between the two shells. 



3.2.1. Implementation 



The Hi and CO datacubes are re-gridded to 2' and 3' spatial grids for the GSH 287+04-17 
and GSH 277+00+36 regions respectively. The fitting is implemented interactively in IDL using 
the non-linear least squares fitting package MPFIT ( |Markwardt ] |2009D , with initial guesses for the 
model function and component parameters either made by eye or based on the fit parameters for 
neighboring spatial points. Examples of fitted spectra are shown in figure [5| 

In cases where the signal component from the shell appears to be a blending of more than one 
Gaussian, the statistical F-test is applied to determine whether the decrease in that arises from 
the addition of the extra component is statistically significant (e.g. Westmoquette et al.||2007 ). We 
require this decrease to be significant at the 4 sigma level for the double fit to be retained. Where 
such double components are present they are labelled according to their velocities, with component 
1 designated as the more blueshifted of the two lines. Post-processing we have become aware of work 
that cautions against this application of the F-test, demonstrating that it is technically statistically 
incorrect despite its extensive use within the community (Protassov et al. 2002). Nevertheless, 



tests on synthetic spectra demonstrate that in our case this approach reliably identifies double 
components, provided that their separation is comparable to or greater than the velocity resolution 
of the data. 



Because we are interested primarily in accurate background subtraction, and because we have 
deliberately chosen regions where the signal emission is well defined, the fitting is free from many 
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of the difficulties that can often plague multi-component Gaussian decomposition (see e.g. Haud 



2000). Although in some cases several solutions may fit the profiles similarly well, or the routine 
may fail to locate the true minimum fit for a given model function, in practice the exact choice of 
background fit does not greatly affect the parameters of the signal component. In the vast majority 
of cases the variation in signal parameters is < 10%, and never greater than ^ 20% even for the 
most ambiguous spectra in our sample. This level of accuracy is sufficient for our needs. 



3.2.2. Column Densities and Number Densities 

Atomic and molecular hydrogen column densities through the shell walls are estimated from 
the fit solutions. TVhi takes values of 1 - 8 x 10^° cm^^, and TVrs of 1 - 60 x 10^° cm-^. At ah 
points where molecular gas is observed, A^H2 ^ ^Hi- 

Here, iVni is estimated from the brightness-column density relation, A^hi = 1.8x 10^^ / 76(HI) dv 



(e.g. Dickey &; Lockman 1990), which holds in the optically thin regime and returns a lower limit 
if optical depth is not negligibly small. A^H2 is estimated using the empirical 'X-factor' relating 
^^CO(J=1-0) integrated intensity and H2 column density, A^H2 ^ ^ J ^KCO) dv, where we adopt 
a value of X = 1.6 x 10^^ K^^ km s"^ cm^^ ( [Hunter et aL]|l997[ ). 



Figure [6] shows Hi column density, A^hi, across the three fitted regions, as well as the ratio 
A^H2/^Hi in regions where molecular gas is detected. Figure [t] shows histograms of A^m, color-coded 
by the presence or absence of spatially coincident CO emission. 

Embedded CO clouds are located in zones of high atomic column density. In contrast, offset 
clouds are associated with very low values of A^hi- This trend is seen clearly in figure [7| where 
the two classes of cloud populate distinct regions within the histograms. It is also reflected in the 
high A^H2/^Hi ratios seen towards offset CO clouds, which peak at ^ 10, ~ 30 and ^ 6 in the 
approaching limb complex, high latitude complex, and molecular drip region respectively. This 
quantitative trend is consistent with the morphological criteria by which the cloud classes were 



originally defined (see ^3.1.3). 



By assuming a depth roughly equal to the characteristic size a feature subtends on the sky, 
a crude estimate of Hi and H2 number densities in the shell walls may be obtained. This yields 
nni ^ 10 cm~^ and ^ 10-100 cm~^ for all three regions. The former is consistent with the 
densities expected of atomic CNM. The latter is somewhat low for molecular clouds, although this 
is likely to be due in part to low fillings factors within the NANTEN beam. 



3.2.3. Velocity Ojfsets Between Atomic and Molecular Emission 



Hi and CO intensity- weighted mean velocities, i;o,Hi and vq^cOi provide information on the 
bulk motions of the atomic and molecular gas. Figure |8] shows images of i^cm, as well as the 
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velocity difference {vq^co — '^o,Hi)- The uncertainty on this velocity difference is small across both 
GSH 287+04-17 regions, where it is estimated to be strictly < 0.2 km s"-*^. However, for GSH 
277+00+36 the lower CO data quality leads to much larger uncertainties of ^ 1 km s~^, and this 
region is therefore ignored in the analysis below. 

No systematic offset is found between the intensity-weighted mean velocities of the two tracers. 
The mean values of {vq^co — '^o,m) averaged over all CO-detected points are consistent with zero for 
both the approaching limb and high latitude complexes. Nevertheless, there is considerable local 
variation on scales of a few parsecs, with (vq^co ~ '^o,Hi) taking values of up to ±2.5 km s~-^. The 
largest of these local velocity differences are observed in embedded CO clouds, and appear to arise 
where CO is either associated with one of an Hi double component pair or offset towards one side 
of a wide atomic line. In contrast, vq^co and i;o,m in offset clouds mostly agree to within ^ ±0.5 
km s~^. For many points this separation is on the borderline of statistical significance. 

With the exception of locally-driven expansion around the Hll region at = (285.9,4.5), 
we are unable to identify any systematic patterns in the Hl-CO velocity offset that might indicate 
coherent macroscopic gas motions such as infall onto individual clouds. 

3.2.4' Velocity Dispersion and Kinetic Temperature 

The velocity dispersion of the Hi gas provides constraints on its kinetic temperature. Under the 
assumption of pure thermal line broadening, Tj^ is determined directly from the Hi Gaussian FWHM, 
Av, meaning that linewidth can be used to obtain a strict upper limit on kinetic temperature, 
Tk.max — 21.9 Ai;^. In the following, Av is defined as the intensity-weighted velocity dispersion, 
cr^, multiplied by the factor y^8/n(2). This corresponds directly to the Gaussian FWHM in the 
case of a single component fit, and gives an equivalent combined weighted velocity dispersion when 
applied to a multiple-component fit. 

Figure [9] shows images of the total FWHM velocity dispersion of the atomic shell emission, 
^'^Hi,fof5 as well as the values for individual Gaussian components, Ai;Hi,i and A'yHi,25 where present. 
Figure [TO] shows histograms of the same quantities, color-coded by the presence or absence of 
spatially coincident CO emission. 

Velocity dispersions are very small; as low as Av^ii^tot ^ 2 km s"-*^ in the narrowest profiles. 
Mean values are 4.4, 4.0 and 4.8 km s~^ in the approaching limb complex, the high latitude complex, 
and the molecular drip region, respectively. Many of the broadest profiles arise from the blending 
of two narrower components. This is illustrated by the fact that replacing the intensity- weighted 
velocity width in blended profiles with Ai;Hi,i or Ai;Hi,2 results in a reduction in these means by 
^ 0.4 km s~^. We also expect there to be more positions where blended profiles remain unresolved 
in our analysis. 

These narrow linewidths indicate cold Hi. For pure thermal broadening, 4 km s~^ corresponds 
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to a kinetic temperature of only 350 K. In practice, turbulence also contributes significantly to 



the observed linewidth, and the true value of Tk is considerably lower. Heiles & Troland (2003) 
find that for the CNM in the solar neighborhood, turbulence and thermal motions generally make 
comparable contributions to Ai;. Their data is best fit by the following expression: log T^^rnax — 
1.32 log — 0.11, where T^, the spin temperature, is equal to the true value of Tj^ for the case of the 
CNM. Applying this relation to our data yields a mean Tj^ of roughly ^ 100 K for all three sample 
sections of shell wall. In non-fitted regions too, there is evidence of similarly narrow linewidths 
throughout both GSH 287+04-17 and GSH 277+00+36. Taken together with the high Hi number 
densities derived above, these linewidths provide robust evidence of shell walls that are dominated 
by cold gas. 

Comparing molecular and non-molecular regions, we find no statistically significant global cor- 
relation between Hi velocity dispersion and the presence of CO. Nevertheless, significant differences 
are observed in Av^.tot towards the two classes of CO cloud. Offset clouds are associated exclu- 
sively with gas with Avm^tot < 4 km s~^, whereas embedded clouds tend to be found towards gas 
with Avm.tot > 4 km s~^. This trend reflects a correlation between the presence of embedded CO 
and multiple blended components in Hi line profiles. 

Finally, CO linewidths are on average narrower than Hi, with mean values of Avco^tot of 
2.0 and 1.9 km s~^. A weak correlation between Av}ii^tot and Avco^tot is found, with Pearson's 
correlation coefficients of ^ 0.5 and ~ 0.3 for the approaching limb and high latitude complexes 
respectively. This trend reflects a preference for wide or blended CO profiles in regions where Hi 
profiles are also wide or blended, and suggests that CO is not strictly localized to single narrow 
zones within large columns of Hi, but exhibits a more complex spatial distribution. 



3.3. Enhanced Molecular Fraction in the Shells 

A critical question for supershell-driven molecular cloud formation is whether we observe an 
enhanced degree of molecularization in the volumes affected by the shells. We here adopt the ratio 
of total CO to total Hi integrated intensity, Yl Ico/ X] ^hi, as a measure of the molecular fraction, 
since both quantities scale approximately linearly with mass. The value of this ratio within the 
volumes affected by our shells is then compared with that in quiescent 'background' regions in 
which no recent supershell activity has taken place. 

Low resolution 15') GASS Hi datacubes are used, since it is necessary to expand the data 
coverage to include regions outside those observed with the ATCA. Because we are concerned with 
integrated intensities summed over large l-h-v volumes the loss of resolution is of no concern. For 



CO, we make use of available NANTEN GPS data (see ^2.2). 



Both shells occupy coherent volumes in l-h-v space, and are seen in all three planes as a 
characteristic void in Hi, often surrounded by a brightened annulus or clumps. We work in the l-v 
plane, defining the extent of the shells by eye at regular latitude intervals throughout the cube. 
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These 2D shell masks are then assembled into a 3D volume enclosing all emission from both the 
shell and the void. All remaining emission in the cube is designated as background. 

There are a number of environmental factors that significantly affect molecular gas abundances; 
most notably altitude above the Galactic midplane, Galactocentric radius, and proximity to spiral 
arms. The spat io- velocity dimensions of each masked cube must therefore be carefully trimmed to 
ensure that the non-shell background regions match as closely as possible the Galactic environment 
in which the shell is found. In practice this involves the following: 1. excluding all emission outside 
the latitude ranges of the shells; 2. excluding velocities containing contributions from local emission; 
3. avoiding the tangent of the Carina Arm (which affects the vicinity around both shells to some 
extent); 4. ensuring that where variation in distance and Galactocentric radius is unavoidable, this 
variation is not only minimal, but also well behaved across the cube and well represented in both 
shell and background regions. The ratio ^ Ico/ ^Hi is then derived for the shell and background 
regions by summing the CO and Hi integrated intensities contained in each. 

The results indicate a significant enhancement in the molecular fraction in both shells. Figure 



11 plots the variation in X^^co/X^^HI with Galactic latitude both in and out of the two shells, 
illustrating that an overabundance of CO is consistently observed throughout almost the entire 
extent of both objects. For GSH 287+04-17 we find a total in-shell ratio of [^Ico/^Im\sh — 
1.3 zb 0.4 X 10"^ compared to a background ratio of /co/ E ^hi]^^ = 0.5 ± 0.3 x 10"^. For GSH 
277+00+36 the equivalent figures are 3.3 + 0.8 x 10~^ and 0.9 + 0.3 x 10~^, implying an enhancement 
in molecular fraction by a factor of ^ 3 for both shells. 

The uncertainties on the above figures are estimated via experimentation with a wide variety 
of different shell masks and spatio- velocity ranges (including some masks defined in the l-h plane). 
A statistically significant increase of E ^CO / E ^hi is consistently measured across all tested con- 
figurations, reflecting a genuine concentration of CO emission in volumes affected by the shells. 
This enhancement spans values of between around 2.0 and 5.5. 

A similar enhancement in mean per-pixel CO intensity is also confirmed over the volumes of 
both shells. While the non-trivial conversion between velocity and distance means that this quantity 
cannot readily be translated into a volume density, it nevertheless demonstrates a genuine increase 
in CO brightness within the two shell regions. This is important, because it demonstrates that the 
enhanced molecular fraction obtained above reflects a genuine over- abundance of CO and does not 
arise from Hi integrated intensity significantly underestimating atomic mass in the CNM-rich gas. 

This result provides clear and direct observational support for supershell-driven molecular 
cloud formation. Crudely, we suggest that as much as 2/3 of the CO emission in both GSH 
287+04-17 and GSH 277+00+36 may arise as a direct result of the influence of these shells. If 
this result is found to hold true on Galactic scales, it will have profound implications for our 
understanding of the role played by stellar feedback on the evolution of the molecular phase. 
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3.4. Total Atomic and Molecular Masses 

It is appropriate to comment briefly on the total swept-up masses of the two shells, since the 
molecular mass of GSH 277+00+36 is newly determined in this work. From the analysis above, the 
total molecular mass of GSH 277+00+36 is estimated to be - 2.1 ± 0.6 x 10^ M©, where the 
error estimate arises from uncertainties on the kinematic distance. The atomic mass is estimated 
by [McClure-Griffiths et ah] (|2000[) as Mri - 3 ± 1 x 10^ Mq. For comparison, [Dawson et"al 



( |2008b[ ) estimate the total Hi and H2 masses of GSH 287+04-17 to be Mri - 7 ± 3 x 10^ M0 and 
^ 2.0 ± 0.6 X 10^ M0, respectively. These fi gures imply a lower degree of molecularization in 
GSH 277+00+36. This is consistent with its location in the outer Galaxy, where ambient densities 
are generally lower. Indeed, the implied initial mean number densities of the pre-shell medium are 
no - 1 cm-i and no - 3 cm-^ for GSH 277+00+36 and GSH 287+04-17, respectively. We note 
that these absolute molecular fractions are distinct from the level of enhancement in the molecular 
fraction in the shell volumes that was computed above. 



4. Discussion 

Both GSH 287+04-17 and GSH 277+00+36 are CNM-rich, with significant masses of associ- 
ated molecular gas. Crucially, the degree of molecularization in the shell volumes appears to be 
higher than in equivalent quiescent regions, strongly arguing for an enhanced level of molecular gas 
production due to the influence of the shells. Yet in an inhomogeneous ISM the ambient medium 
is pre-structured; containing zones of cold, dense material whose origin predates the supershells. 
Molecular gas may be formed in-situ from the swept-up medium, but this process must occur in 
tandem with the interaction of the shells with pre-existing molecular clouds. 

The embedded and offset CO clouds identified in §3.1| provide some of the best candidates for 
these two processes. Embedded clouds - with their close association with the material of the main 
shell walls - are strongly suggestive of molecular gas that has condensed from the swept-up atomic 
medium. In contrast, offset clouds - located at the tips of swept-back Hi features - may be well 
explained as the remnants of an initial medium that was peppered with dense clumps. 

The remainder of this paper demonstrates the plausibility of these two scenarios for archetypal 
examples of embedded and offset clouds. We use these objects as a basis from which to explore some 
key issues surrounding molecular cloud formation and destruction, focussing particularly on cloud 
sweep-up/formation timescales, the outcome of shock-cloud interactions, and the requirements for 
shielding of the CO molecule. 

In reality there is unlikely to be a pure dichotomy between the in-situ formation of new molec- 
ular gas from a purely ambient atomic medium and the interaction of the shell with a population of 
fully pre-existing molecular clouds. Nevertheless, this idealization provides a meaningful starting 
point from which to explore the relationship between supershells and the molecular ISM. 
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Finally, the evolution of the molecular ISM in supershells links directly to the question of 
the dominant processes driving star formation within the Galaxy. We round off our discussion 
by making some brief remarks on star-formation activity within the shell clouds, prior to further 
dedicated study on this topic. 



4.1. In-Situ Molecular Cloud Formation: Embedded CO Clouds 



Supershells enable the rapid production of molecular gas by increasing both the number density 
and column density of the atomic medium in localized regions. The strong dependence of molecule 
formation rates on density means that the formation of molecules should proceed rapidly in swept- 
up shells relative to the background, low-density medium. Furthermore, the large spatial scales 
over which supershells operate, and the sheer quantities of material they are able to accumulate 
over their lifetimes, lead to high column densities within shell walls. This provides shielding from 
the background UV field, which is essential to the survival of molecular species. 

The double filament structure observed in the approaching limb complex of GSH 287+04-17 
(figure [2] panels al and a2; 285.3 < / < 286.4°, 4.2 < 6 < 6.0°) is a particularly striking example 
of a feature well explained by the in-situ formation of molecular gas in the swept-up shell. The 
embedded CO clouds form an integral and coherent part of the curved Hi wall, with both tracers 
closely associated in l-b-v space. High Hi column densities suggest reservoirs of raw material 
available for molecule formation, as well as for shielding of the CO molecule. 

Numerical studies suggest that a minimum (ID) extinction of Ay ^ 0.7 is required for the 



rapid formation and sustained survival of molecular gas from the atomic medium (Bergin et al. 
2004| [Glover et aL]|2010[ [Wolfire et aL]|2010|). The relation Ay = A^hi/(1.9 x 10^^ cm^^) (e.g. 



Draine &; Bertoldi 1996) may be used in conjunction with the Hi column density measurements 



from §3.2.2| to estimate Ay in the material surrounding the CO clouds. This returns values of 
(Ay) ^ 0.2 towards lines of sight coincident with embedded CO. However, since A^hi is a lower 
limit in the case where optical depth is not negligible, this figure is expected to miss a significant 
fraction of the CNM that dominates the shell walls. Similarly, it cannot account for CO-dark H2, 
which requires much lower extinctions in order to form, and which must exist in some quantity in 
envelopes around CO-bright clouds. The contribution from both of these dark fractions may be 
considerable, and can range from several 10% to several times the CO-bright molecular column 
density ( [Grenier et al^|2005l [Wolfire et al.||20Tol ). 



Infra-red excess measurements provide a means of evaluating this 'dark gas' fraction. Under 
the assumption of a constant Galactic gas-to-dust ratio, the dust opacity at wavelength A, r^, is 
directly proportional to the total hydrogen column density along the line of sight. We make use of 



reprocessed IRAS 60 jam and 100 jj^m images from the IRIS project (Miville-Deschenes &; Lagache 



2005) to evaluate the 100 fim opacity, rioo- This is then compared with our Hi and CO datacubes to 



obtain linear relationships between rioo and Hi and CO integrated intensities, and hence to assess 
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the excess column density contribution arising from material not accounted for by either of these 
two transitions. 

We work with an 8° x 6° region centered on the GSH 287+04-17 approaching limb complex, 



and follow the prescription of Douglas &; Taylor (2007) in order to obtain a measurement of the 



excess opacity Tioo,ex- The present analysis differs from their treatment in the following points: 1. 
We do not explicitly account for the contribution from ionized material, which is generally negligible 
across the region. 2. We do not apply their 'correction factor' that modifies A^hi to include an 
rough estimate of the contribution from cool gas, since it is explicitly this contribution that we 
wish to assess. 



The results of this analysis are shown in Figure [T2j Significant excess opacity is observed in the 
entire section of shell wall, concentrated especially in regions of embedded CO. The highest values 
of noo,ea: occur in the vicinity of the small Hll region at (/, b) = (285.9°, 4.5°). These are assumed 
to include a significant contribution from ionized material, and are therefore ignored. However, 
even outside of this anomalous region, Tioo,ea: takes consistently high values, averaging ^ 9 x 10~^ 
for CO-detected positions. 

The column density implied by this value is estimated by applying the scaling factor relating 
Tioo and Hi column density that was obtained during the above analysis. This factor is rioo = 
(4.7 ±0.2) X 10"^^ Nm, which implies missing hydrogen column densities of ^ 2 x 10 cm . This 
in turn translates to Ay ^ 1, which combined with the contribution directly traced in Hi, gives a 
total extinction of ^ 1.2. Strictly speaking, this figure applies to all material along the line of sight 
outside of the CO-bright core region, implying a ID extinction at the surface of the CO cloud of 
0.6. However, the parsec-scale resolution of the datasets ensures that this figure underestimates 
the true peak values of Ay in the region. These results imply that enough material exists in the 
swept-up wall to permit the existence of a stable core of CO-rich molecular gas. 

The mass of the entire complex is estimated as ^ 2.5 x 10^ Mq. Of this, less than half is 
traced by Hi and CO, and the remainder is inferred from IR excess, demonstrating the importance 
of material invisible to these two standard tracers. By assuming that the section of wall possesses a 
depth somewhere between the minimum and maximum widths it subtends on the sky, an estimate 
is obtained of the average initial number density required to form the complex via the sweeping 
up of the pre-shell medium. Taking the center of expansion of the shell as (/,6) = (287.5°, 3.0°) 
( Dawson et al.|[2Q08b ), results in an initial number density of (nn) ^3 — 10 cm~^. This is slightly 



enhanced with respect to the mean density of the pre-shell medium calculated over the entire shell 



volume from Hi and CO alone (see ^3.4). However, it is consistent with a scenario in which initial 
over-densities within the ambient medium play a role in determining the sites at which efficient 
molecular cloud formation can occur. 

Finally, it is essential to consider the the formation of molecular gas in supershells in the 
context of current cloud formation theory. Modern theory favors a picture in which clumpy sheets 
and filaments of CNM are formed through the compression, cooling and fragmentation of WNM 
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in large-scale shocks and colliding flows, of which the case of a supershell expanding into the 



surrounding medium is but one example (e.g. Hennebelle &; Perault 1999; Audit & Hennebelle 
2005; Vazquez-Semadeni et al.||2007 ; Heitsch fc Hartmann||2008 ). Densities in these fragments may 
greatly exceed canonical values for the CNM, as a consequence both of the global ram-pressure and 
the presence of transient turbulence-driven density enhancements. This increased density leads to 
fast chemical timescales, and provided that the interaction zone is able to accumulate sufficient 
material, chemical models suggest that the transition to molecular gas may occur on timescales 
of - 10^"'^ yr dKoyama fc Inutsuka|[2000{ [Bergin et aT]|2004| [Heitsch fc Hartmann|[2008[) . In this 



context it is interesting to note the correlation of embedded CO with complex Hi velocity profiles 



(see ^3.2.4), which may indicate the stirring up of material in interaction zones, or the presence of 



further local compressive motions within the shell wall, either of which could act to further increase 
density in localized regions. 

For the conditions appropriate to our objects, the colliding flows model is able to form dense 
clouds on considerably shorter timescales than alternative methods such as the purely gravitational 



fragmentation of the swept up shell (e.g. Elmegreen et al. 2002|). ID models of molecular cloud 
formation that explicitly include the chemistry of CO and H2 formation suggest that a flow velocity 
of as little as ^ 10 km s~^ can produce large quantities of CO-rich molecular gas in ~ 10^ years 



for an initial number density of a couple of atoms per cubic centimeter ( Bergin et al.||2004 ). These 
figures are consistent with the estimated age and expansion velocity of GSH 287+04-17 as a whole, 
and formation timescales should be even shorter for the slightly enhanced initial densities implied 
for the specific section of shell wall discussed above. For GSH 277+00+36, which is both older 
and more quickly expanding, the conditions for formation are also met easily. Moreover, these 
formation timescales are expected to become even lower when models are expanded to include 
reahstic, turbulent gas dynamics (Glover &; Mac Low 2007; Glover et al.||2010 ). 



Magnetic fields are expected to play an important role in this cloud formation process (Inoue 



& Inutsuka 2008, 2009; Heitsch et al. 2009). Magnetic support acts to oppose the formation of 



very dense gas, leading to the additional requirement that the component of B perpendicular to 
the flow direction must be vanishingly small in order for molecule formation to proceed. This 
rather stringent requirement may be crucial in pre-selecting the sites where molecular clouds can 
form, and is pleasingly consistent with the observational reality that supershells contain at best a 
scattering of molecular clouds. Although at present we have no information on the magnetic fields 
in either of our shells, this subject should motivate further observational study. 



4.2. Offset CO Clouds: Pre-Existing Molecular Gas 

The approaching limb complex of GSH 287+04-17 also contains several promising candidates 
for pre-existing molecular clouds. Between 285.5 < / < 288.5°, 5.4 < b < 6.2° (figure [2] panel ai), 
CO-tipped Hi 'fingers' point radially inward toward the shell interior, with the main walls pulled 
back, arch-like, around them. 
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In contrast to the embedded clouds described above, these smah offset CO-clouds are preferen- 



tiahy associated with very low values of A^hi (see ^3.2.2). Similarly, there is no significant IR excess 



coincident with the clouds (see figure 12), indicating that there is little or no dark gas undetected 
in Hi or CO. Correspondingly, not only is there minimal material for shielding, but there is also no 
substantial reservoir of raw material available for molecular gas formation. The atomic to molec- 
ular mass ratio in CO-detected pixels averages only ^0.1, and the mean extinction arising from 
material not traced by CO is estimated from A^hi to be Ay ^0.1. This is an order of magnitude 
smaller than required to shield CO. We conclude that the ongoing formation of molecular gas in 
these clouds in their present configuration is not viable. 

A simple check may be performed on a scenario in which the present CO clouds are remnants of 
molecular gas that pre-dates the shell, by comparing theoretical cloud destruction timescales with 
the survival times implied by the data. A lower limit to the cloud survival time may be estimated 
from the shell expansion velocity (^ 10 km s~^) and cloud-to-main- wall distances (20 — 40 pc), 
by assuming that the passage of the shell has not accelerated the clouds significantly and ignoring 
projection effects. This implies that the interaction between the clouds and the wall occurred a 
minimum of 2 — 4 Myr ago. Similar timescales are obtained for the offset clouds in the molecular 
drip region of GSH 277+00+36 (Figure |4) panels al and a2). Here, cloud-to -main-wall distances 
of 20 — 30 pc and an expansion velocity of ^ 20 km s~^ imply survival times of at least 1 — 1.5 Myr. 

Dense clouds interacting with an expanding shell will be subject to dynamical disruption from 
the shell-cloud interaction, followed by thermal evaporation by the hot interior medium in which 
they become entrained. In addition, the physical stripping and dissipation of cloud material also 
renders them increasingly vulnerable to the UV dissociation of CO. The observable lifetime of the 
entity we see as a CO cloud may therefore be shorter than the survival time of the dense material 
itself. 

The dynamical interaction of a dense cloud with a shocked flow is commonly parameterized 
in terms of the cloud crushing time, tec = (^o/'^i)(^d/^i)^'^5 where tq is the cloud radius, Vi is the 
velocity of the shock in the ambient medium, and rid ^iid Ui are the number densities of the cloud 



and the ambient medium respectively (e.g. Klein et al. 1994). For the case of a section of shell 
wall impacting a pre-existing molecular cloud, we adopt the parameters rg = 2 pc, Vi = 10 km s~^, 
rici = 100 cm~^ and rii = 10 cm~^, which are appropriate for the shell and clouds in their present 
state, and correspond to tec ^ 0.6 Myr. For the ideal adiabatic case clouds are typically destroyed 
on timescales of a few to ^ 10 tec ( Klein et al.|1994 ; Nakamura et al.|2006 ). However, for the present 
case radiative cooling cannot be ignored over the timescales of the interaction. Radiative cooling 
tends to inhibit cloud destruction, encouraging the formation of over-dense clumps and filaments. 



and potentially extending survival times significantly (Mellema et al. 2002; Orlando et al. 2005) 



This suggests that dynamical disruption of pre-existing molecular clouds will not be so fast as 
render them unobservable on the relevant timescale of 1 — 4 Myr 



For a shell thickness of ^ 10 pc, the cloud will pass through the shell and into the hot interior 
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on timescales of 1 to several Myr, depending on the drag efficiency. In the interior regime the 
material flowing past it will be more diffuse and hence less dynamically disruptive, but thermal 
conduction may become important. The classical evaporation tim escale for clouds embe dded in 
fully ionized medium is given by teyap ^ 3.3 x 10^^ rid T- ^^^r^^ yr (Cowie & McKee 1977), which 
equates to ^ 10^ yr for the present clouds, assuming Ti = 10^ K. Thermal evaporation is therefore 
unlikely to be important over the lifetime of the shell. 

The CO destruction rate by UV may be expressed as 

I'coG'ofco e-^^ooAv s-\ where /qq is the 
unshielded CO photodissociation rate in the Draine field, Gq is the ratio of the incident radiation 
field to the Draine field, fco is the CO self shielding factor, and bco is the dust extinction coefficient 
for CO (see Appendix C of Wolfire et al.||2010 ). Values of Iqq and fco are taken from Visser et al. 
( 2009| ), where fco is estimated from their table 7, using values of A^H2 from our data and assuming 
a standard H2 to CO abundance ratio of 10^. Following Wolfire et al. (2010), bco is taken as 3.2, 



Gq is assumed to be 1, and Ay is estimated as 0.8 based on the values of A^H2 ^Hi at the 
spatial centers of the CO clouds. This gives an estimated hfetime of rco ^ 5 Myr for CO at the 
center of the clouds in their current form, implying that at least some molecular material should 
remain observable for several million years. In contrast, the survival time at the surface of the CO 
clouds where Ay ^ 0.1 is as little as ^ 10^ yr. This is consistent with a scenario in which the 
dynamical stripping of material from around an originally stable CO cloud renders it increasingly 
vulnerable to UV dissociation. 

While all of the timescales estimated above are necessarily rough, they nevertheless illustrate 
the plausibility of a scenario in which many of the offset clouds we observe in our two shells are 
remnants of pre-existing molecular material. However, it must be noted that not all offset clouds in 
the present sample are so neatly explained by this simple scenario. Particularly those at very high 
latitudes, for which it is difficult to postulate the a-priori existence of a large quantity of pre-existing 
molecular gas. A more realistic description of supershell evolution might contain 'hybrid' scenarios, 
in which clumpy shell walls containing newly-formed CO clouds are carved by later episodes of 
energy input into the configurations now observed. 

Finally, we note that we have not commented on the possible role of Rayleigh- Taylor instabil- 
ities in creating shell structure, which has been suggested as a possible means of forming the Hi 



'drips' observed in GSH 277+00+36 (McClure- Griffiths et al. 2003). The presence of pre-existing 



dense structure in the inhomogeneous ambient medium is likely to affect the development of the 



instability, modifying the properties of the RT fingers that form (e.g. Jun et al. 1996). However, 
a full consideration of the interplay between instabilities in the shell wall and pre-existing dense 
clouds is beyond the scope of this paper. 
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4.3. Remarks on Star Formation 



Since molecular clouds are the sites of star formation, the creation and destruction of molecular 
material in supershells is of potentially critical importance to Galactic star formation rates. 



Fukui et "al] ( 1999| ) compared the molecular cloud population of GSH 287+04-17 with young 



stellar object (YSO) candidates extracted from the IRAS point source catalogue, finding a number 
of convincing correlations that strongly argued for ongoing star formation in the shell clouds. 



Figure [T3| reproduces this source list for GSH 287+04-17, while Figure [T4| extends the analysis to 
GSH 277+00+36. The present high-resolution comparison of CO and Hi allows us to place these 
potential star forming regions in their proper context within the ISM of the shell walls. 

The molecular clouds of the GSH 287+04-17 approaching limb complex show excellent cor- 
relation with a number of IRAS sources. Particularly striking is the presence of YSO candidates 
in all three of the small offset CO clouds discussed in §4.2[ As promising examples of molecular 
material pre-dating the shell, this raises the possibility that star formation is being triggered in 
these clouds as a direct result of their interaction with the shell. Indeed, in shock-cloud interactions 
where radiative cooling is important, compression and cooling may lead to localized gravitational 
collapse even as much of the mass of the cloud is stripped away ( [Foster fc Bossl|1996t [Fragile et aL| 
2004; van Loo et al.||2007 ). Similarly, the larger of the two offset clouds in the molecular drip region 
of GSH 277+00+36 is also associated with a YSO candidate. 

Perhaps more importantly, the embedded molecular cloud in the approaching limb complex 



of GSH 287+06-17 discussed in ^4.1 contains an active massive star forming region at 



(285.9,4.5). This region contains a small optical Hll region (see Figure 12), two luminous IRAS 



YSOs, and the recently confirmed 2MASS cluster DBSB 49 ( [Dutra et all2003 [). The most massive 
member of this small open cluster is of spectral type BOV, and it has an estimated age of 2.1 ± 0.3 
Myr ( Soares et al.|2008 ), placing the onset of star formation relatively recently in the shell's history. 
Given this cloud's proposed origin as an object formed in-situ in the shell wall, this suggests the 
triggered formation of both molecular gas and stars in ISM flows on total timescales of ^ 10^ 
yr. This is highly consistent with a picture in which the onset of star formation in locally self- 
gravitating clumps occurs rapidly after sufficiently dense conditions are reached for the creation of 
the parent molecular cloud ( Hartmann et al.|[2001 ; Vazquez-Semadeni et al. 2007). 



Embedded molecular material in GSH 277+00+36 shows similarly strong associations with 
luminous IRAS sources, as well as other clear indicators of star formation activity. In particular, 
the giant Hll region RCW42 (/=274.0°, 6=-l.l°) is located within the prominent embedded CO 
complex discussed in §3.1.2[ attesting to active massive star formation in the shell walls. This 
region also contains the deeply embedded 2MASS stellar cluster DBSB 38 ( Dutra et al.||2003 |, and 
several other luminous YSO candidates associated with the wider CO cloud complex. In addition, 
several other close associations between IRAS sources and CO clouds in both shells are shown in 



Figures 13 and 14 
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Taken together, these results attest to the ongoing formation of stars throughout the molecular 
material in both GSH 287+04-17 and GSH 277+00+36. We suggest that future studies might 
build on these preliminary findings in order to better quantify the global star formation rates in 
supershell-associated molecular clouds, and explore their wider role in the Galactic ecosystem. 

5. Conclusions 

We have presented parsec-scale resolution observations of Hi and ^^C0( J=l-0) in two Galactic 
supershells, GSH 287+04-17 and GSH 277+00+36, with the aim of investigating the role played 
by supershells in the evolution of the molecular ISM. The main findings may be summarized as 
follows. 

1. Both shells contain large quantities of associated molecular gas in the form of discrete co- 
moving CO clouds distributed throughout the atomic shell walls. These high resolution observations 
reveal rich substructure in both tracers. Molecular gas is seen elongated along the inner edges of 
atomic shell walls, embedded within Hi filaments and clouds, or taking the form of small CO clouds 
at the tips of tapering 'fingers' of Hi. We note a similarity in features observed in both objects, 
despite differences in location and evolutionary stage. 

2. The atomic shell walls are dominated by cold gas, showing narrow linewidths reaching as 
low as Av ~ 2 km s"-"^. Mean temperatures and densities are estimated to be roughly ^ 100 K 
and no ~ 10 cm~^ in regions where the spectral signature from the shells is well determined. 

3. An enhanced level of molecularization is observed over the volumes of both GSH 287+04-17 
and GSH 277+00+36, providing the first direct observational evidence of increased molecular cloud 
production due to the influence of supershells. Our results imply that the amount of molecular 
matter in the volumes affected by these shells is enhanced by ^ 3 times with respect to neighboring 
regions. If this is found to hold true on Galactic scales, it has powerful implications for our under- 
standing of the role played by stellar feedback on the evolution of the molecular phase. Already, the 
confirmation of this phenomenon in two shells at quite different locations and evolutionary stages 
is very compelling, and we suggest that more foUowup work should be undertaken to repeat this 
analysis for other objects. 

4. CO clouds embedded in the main atomic shell walls provide excellent candidates for the 
in-situ formation of molecular gas from the swept up medium. This scenario has been explored 
in detail for archetypal examples of embedded molecular clouds, demonstrating that the formation 
timescales implied by theory are consistent with shell ages and number densities, and that the 
requirements for the shielding of the CO molecule are met. We thus confirm on a local scale the 
viability of the triggered formation of molecular clouds due to the influence of the shells. 

5. Small offset CO clouds located at the tips of tapering 'fingers' of Hi may be the remnants 
of molecular gas present in the ISM prior to the formation of the shells. We have demonstrated the 
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plausibility of this scenario for archetypal examples of such offset clouds, by comparing estimates 
of cloud destruction times with the survival times implied by the data. 

6. A preliminary examination of YSO candidates in the shell regions confirms that active 
star formation is occurring in the molecular clouds in both GSH 287+00-17 and GSH 277+00+36. 
This includes massive star forming regions in embedded molecular cloud complexes, as well as less 
luminous YSO candidates associated with offset CO clouds. 
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Velocity; -33,46 km/s Velocity: -30.99 km/s 




Galactic Longitude Galactic Longitude 

Fig. 1. — Channel maps of the GSH 287+04-17 region. Each panel shows mean brightness tem- 
perature over a velocity interval of 2.47 km s~^. Greyscale images are Hi data. Pink contours show 
^^CO(J=1-0), beginning at a starting level of 0.7 K and incremented in steps of 1.0 K. Dotted lines 
mark the boundary of the area observed in ^^CO(J=1-0). The white dashed line in panel g marks 
the approximate outline of the main body of the shell. 
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Fig. 1. — cont. 
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Fig. 2. — Integrated intensity images showing subregions of interest in GSH 287+04-17. Greyscale 
images are Hi and pink contours are ^^CO(J=1-0). The velocity integration ranges and contour 
levels are as follows: (al) —26.5 < visr < —19.9 km s"-*^, 1.5+3.0 K km s"-*^; (a2) —23.2 < visr < 
-21.5 km s-\ 1.3+1.0 K km s-^ (b) -33.0 < visr < -25.6 km s-\ 1.5+5.0 K km s-^; (c) 
— 14.1 < Visr < —10.8 km s"-*^, 1.5+3.0 K km s"-*^. Panels (al) and (b) show the regions referred to 
in the text as the 'approaching limb complex' and 'high latitude complex', respectively. CO clouds 
labelled 1 and 2 indicate those specifically referred to as 'embedded' and 'offset' in the text. 
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Galactic Longitude Galactic Longitude 

Fig. 3. — Channel maps of the GSH 277+00+36 region. Each panel shows mean brightness tem- 
perature over a velocity interval of 2.47 km s~^. Greyscale images are Hi data. Pink contours show 
^^CO(J=1-0), beginning at a starting level of 0.8 K and incremented in steps of 1.0 K. Dotted lines 
mark the boundary of the area observed in ^^CO(J=1-0). 
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Fig. 3. — cont. 
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Fig. 3. — cont. 
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Fig. 4. — Integrated intensity images showing subregions of interest in GSH 277+00+36. Greyscale 
images are Hi and pink contours are ^^CO(J=1-0). The velocity integration ranges and contour 
levels are as follows: (al) 40.8 < visr < 43.3 km s~\ 1.5+1.0 K km s~^ (a2) 42.5 < visr < 4:4:.! km 
s-^ 1.2+1.0 K km s-^ (b) 40.8 < visr < 43.3 km s-\ 1.5+2.0 K km s"^; (c) 43.3 < visr < 48.2 
km s~^, 1.5+1.0 K km s~^. Panel (al) shows the region referred to in the text as the 'molecular 
drip region'. The spectrum shows CO emission from the single-pixel cloud in panel a2. CO clouds 
labelled 1 and 2 indicate those specifically referred to as 'embedded' and 'offset', in the text. 
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Fig. 5. — Examples of Gaussian decomposition of Hi and ^^CO(J=1-0) spectra in the GSH 287+04- 
17 region. Solid black lines show observed spectra. Dashed black lines show Gaussian components 
fitted to signal emission arising from the shell walls. Dashed grey lines show components fitted to 
non-related background emission. 
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Fig. 6. — A^Hi (panels labelled 1) and A^h2/^hi (panels labelled 2) derived from Gaussian fits 
to spectra in selected regions of the two supershells. Panels a and b show the approaching limb 
and high latitude complexes of GSH 287+04-17. Panel c shows the molecular drip region of GSH 
277+00+36. Black contours are A^h2 beginning at a starting level of 2.4 x 10^^ cm~^ and incremented 
in steps of 8.0 x lO^^ cm-^. 
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Fig. 7. — Histograms of A^hi derived from Gaussian fits to spectra in selected regions of the 
two supershells. The charts are color-coded according to the presence or absence of ^^CO(J=1-0) 
emission at a given position. Panels a and b show the approaching limb and high latitude complexes 
of GSH 287+04-17. Panel c shows the molecular drip region of GSH 277+00+36. 
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Fig. 8. — Hi intensity- weighted mean velocity, 'Uo,hi (panels labelled 1), and the velocity offset 
between Hi and -'^^CO(J=1-0), ('^0,00 — '^o,hi) (panels labelled 2), derived from Gaussian fits to 
spectra in selected regions of the two supershells. Maps a and b show the approaching limb and high 
latitude complexes of GSH 287+04-17. Map c shows the molecular drip region of GSH 277+00+36. 
Black contours are A^h2 beginning at a starting level of 2.4 x 10^^ cm~^ and incremented in steps 
of 8.0 X 10^0 cm-2. 
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Fig. 9. — Hi velocity dispersion, derived from Gaussian fits to spectra in selected regions of the 
two supershells. Panels a and b show the approaching limb and high latitude complexes of GSH 
287+04-17. Panel c shows the molecular drip region of GSH 277+00+36. The main maps show 
the quantity A'UHi^tot — \/8/n(2) ay^toti where ciy^tot is the intensity- weighted standard deviation of 
all shell wall emission at a given position. For positions where the shell emission was best fitted by 
two Gaussian components, inset panels show the velocity dispersions, Avm,! and A^'hi,2 of these 
components separately. These locations are also marked with white crosses. Contours are A^h2 5 
starting at a level of 2.4 x 10^^ cm~^ and incremented in steps of 8.0 x 10^^ cm~^. 
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Fig. 10. — Histograms of Av^^i derived from Gaussian fits to spectra in selected regions of the two 
supershells. Grey shaded charts show Avm^tot^ color-coded according to the presence or absence of 
^^CO(J=1-0) emission. Solid, dashed and dotted lines delineate the equivalent regions for Avm,!- 
The top axis shows the upper limit on kinetic temperature, given by Tk^max — 21.9 Av'^. Panels a 
and b show the approaching limb and high latitude complexes of GSH 287+04-17. Panel c shows 
the molecular drip region of GSH 277+00+36. 
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Fig. 11. — Variation of [J2 Ico/J2 ^ailsh and Ico/J2 ^ailbg with Galactic latitude. The spatio- 
velocity ranges of the sub-cubes used in this analysis are as follows: GSH 287+04-17: 284.0 < 
/ < 293.0°, 0.83 <b < 8.17°, -40.0 < visr < -10.0 km s-^; GSH 277+00+36: 270.0 < / < 282°, 
-2.27 < b < 1.47°, 20.0 < visr < 50.0 km s-^. Plotted [Elco/Elm]sh and [Elco/Elm]bg 
data points are calculated by summing the total CO and Hi integrated intensities over the defined 
shell and background regions within slices each spanning 0.33° in latitude. 
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Fig. 12. — Top panel: The greyscale image shows the excess dust opacity at 100 /xm, corresponding 
to gas not traced in either Hi or CO. Black contours are drawn at the 5a detection level of 2.5 x 10~^, 
and incremented every 2.5 X 10"^ to a level of 1.5 x lO""^. White contours plot SHASSA Ha emission 



(Gaustad et al. 2001) at a starting level of 550 dR and incremented in steps of 100 dR. Bottom 



panel: Greyscale Hi image and CO contours as in panel al of figure [2j 



- 41 - 




Fig. 13.— IRAS YSO candidates in the GSH 287+04-17 region. The main panel shows ^^C0{J=1- 
0) integrated over the entire sheh velocity range (pink contours) and IRAS YSO candidates (black 
crosses). Contour levels are 1.5 + 5.0 K km s"-*^. Panels (a) and (b) show greyscale Hi integrated 
intensity in selected subregions overlaid with CO contours at levels of 1.5 + 3.0 K km s~^. Velocity 
integration ranges are (a) —26.5 < visr < —19.9 km s"-*^ and (b) —14.1 < visr < —10.8 km s"-*^. 
IRAS luminosities are calculated only for those objects coincident with molecular gas, for which a 
genuine association with the shell may therefore be assumed. The IR luminosity in the four IRAS 



bands is converted to an estimated bolometric luminosity using the correction factor of Myers et al. 



(1987), under the assumption of an SED that peaks at 100 jj^m. Filled crosses show those sources 
for which the la uncertainty on L^^/ is < 20%. Dotted crosses show those for which the uncertainty 
is greater, or for which the derived luminosity is an upper limit. 
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Fig. 14.— IRAS YSO candidates in the GSH 277+00+36 region. The main panel shows 12C0(J=1- 
0) integrated over the entire sheh velocity range (pink contours) and IRAS YSO candidates (black 
crosses). Contour levels are 2.0 + 5.0 K km s"-*^. Panels (a), (b) and (c) show greyscale Hi integrated 
intensity in selected subregions overlaid with CO contours. Velocity integration ranges and contour 
levels are (a) 40.8 < visr < 43.3 km s-\ 4.0+5.0 K km s-\ (b) 29.3 < visr < 36.7 km s-\ 1.5+1.0 
K km s"-*^, and (c) 36.7 < visr < 43.3 km s~-^, 3.0+5.0 K km s~-^. IRAS luminosities are calculated 
only for those objects coincident with molecular gas, for which a genuine association with the shell 
may therefore be assumed. The IR luminosity in the four IRAS bands is converted to an estimated 



bolometric luminosity using the correction factor of Myers et al. (1987), under the assumption of 



an SED that peaks at 100 fim. Filled crosses show those sources for which the la uncertainty on 
Li,ol is < 20%. Dotted crosses show those for which the uncertainty is greater, or for which the 
derived luminosity is an upper limit. 



